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Plasmonic crystals for ultrafast nanophotonics: 
Optical switching of surface plasmon polaritons 
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We demonstrate that the dispersion of surface plasmon polaritons in a periodically perforated gold 
film can be efficiently manipulated by femtosecond laser pulses with the wavelengths far from the 
intrinsic resonances of gold. Using a time- and frequency- resolved pump-probe technique we observe 
shifting of the plasmon polariton resonances with response times from 200 to 800 fs depending on 
the probe photon energy, through which we obtain comprehensive insight into the electron dynamics 
in gold. We show that Wood anomalies in the optical spectra provide pronounced resonances in 
differential transmission and reflection with magnitudes up to 3 % for moderate pump fluences of 
0.5 mj/cm 2 . 

PACS numbers: 73.20. Mf, 78.47.J-, 78.66.Bz, 78.67.Pt 



Nowadays plasmonics attracts much research inter- 
est in nanophotonics inspiring scientists to develop a 
new paradigm in data processing based on nanometal- 
lic circuitry 1,2 . The key object of plasmonics is a surface 
plasmon polariton (SPP) - a coupled oscillation of the 
electromagnetic field and the electron plasma in metals 3 . 
Excitation of a SPP leads to significant electromagnetic 
energy localization near the metal-dielectric interface, 
thereby enhancing nonlinear effects and light-matter in- 
teraction. Current state-of-the-art in telecommunica- 
tions requires plasmonics to be active, i.e. a possibil- 
ity for control by means of an external stimulus must be 
provided on the order of a few nanoseconds or shorter 4 . 

One of the approaches satisfying the strict criteria is 
magnetoplasmonics employing the influence of an exter- 
nal magnetic field on the SPP propagation constant 5,6 . 
While the modulation efficiency in magnetoplasmonic 
structures may be as large as tens of percent the oper- 
ation rate is limited by the magnetization dynamics oc- 
curring on the nanoseconds timescalc. Transient changes 
in the optical properties of plasmonic structures can be 
also achieved via application of intense femtosecond laser 
pulses 4 ' 7-14 . Here, the real and imaginary parts of the di- 
electric constant of a metal change in the matter of sev- 
eral hundreds of femtoseconds 15-17 . Since the propaga- 
tion constant of SPPs is determined by the permittivity 
of metal and dielectric this opens new horizons for ultra- 
fast control of SPPs. However, the high reflectivity of 
smooth gold surfaces prevents most of the incident elec- 
tromagnetic energy from absorption in the metal. Here, 
SPPs can be exploited to provide electromagnetic energy 
localization near the metal-dielectric interface and con- 
sequently to increase the energy absorption in gold 18 ' 19 . 

Recently the optical response of plasmonic crystals 
with periodically perforated gold or perforated dielectric 
on top of gold was investigated using a two-color pump- 
probe technique with the pump photon energy far from 



the SPP resonances 11,12 . It was demonstrated that the 
differential reflectivity of the probe beam tuned in res- 
onance with a SPP can be modulated up to 10% on a 
sub-picosecond timescale. Large values of the modula- 
tion parameter were accomplished by using intense laser 
pulses with fluence $ ~ 50mJ/cm 2 and by tuning the 
resonances of the probed SPPs close to the d-band tran- 
sitions in gold (hujii, = 2.3 cV 16 ). However, the drawback 
of this approach is the large absorption in the vicinity of 
the d-band transition resulting in a huge intrinsic SPP 
damping. A similar obstacle, but with even larger opti- 
cal losses, is present in Al plasmonic crystals 20 . 

Here, we demonstrate that the dielectric constant of 
the metal in a plasmonic crystal can be efficiently modi- 
fied at femtosecond time scale via SPP excitation, leading 
to ultrafast modification of the SPP dispersion. This is 
possible in the spectral range far from the intrinsic res- 
onance frequencies of the metal system. From the tem- 
poral and spectral dependence of differential reflectivity 
and transmission we extract comprehensive information 
on the electron energy relaxation in gold. We achieve 3 % 
probe intensity modulation with about an order of magni- 
tude smaller pump fluence of $ w 0.5 mJ/cm 2 compared 
to previous works 12 . 

The investigated plasmonic crystal is schematically 
shown in Fig. la. It comprises a periodically perforated, 
120-nm-thick gold film deposited on top of a smooth 
ferromagnetic layer of rare-earth iron-garnet BUG with 
2.5 fim thickness. The period is d = 595 nm with a groove 
depth of 120 nm and an air groove width of 110 nm [see 
A]. The time- resolved differential transmission and re- 
flection were measured by a pump-probe technique using 
a Ti-Sa laser with a pulse repetition rate of about 80 
MHz, center wavelength around 800 nm and a pulse du- 
ration of tjj = 30 fs [see B] . In the experiment the pump 
beam entered the sample in the plane perpendicular to 
the slits, while the probe beam was incident in the slit 
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FIG. 1: (a) Scheme of the plasmonic crystal for active con- 
trol of SPPs. (b) Dispersion of the SPP propagating at the 
BilG/Au [(1),(2) dashed curves] and at the Au/air [(3) dashed 
curve] interfaces calculated by the S-matrix method. Black 
solid lines correspond to the dispersion of the plane waves 
in free space incident at 0°, 12°, 17°, and 22° angle. (c,d) 
Experimentally measured transmission (solid lines) and re- 
flection (dashed lines) spectra for the pump (c) and probe (d) 
configurations. Angle of incidence is 17°. Incident light is 
p-polarized in (c) and s-polarized in (d). 



plane (Fig. la). The incidence angles 9\ and #2 for both 
pump and probe were set to 17° with an aperture angle 
around 6°. The plasmonic crystal was designed to have 
SPP resonances in this configuration in the energy range 
1.5-1.7cV for both the pump and probe beams. 

The dispersion diagram for the SPP modes in this plas- 
monic crystal calculated with the scattering matrix (S- 
matrix) method 21 ' 22 reveals that the pump beam excites 
the SPPs in the 2-nd band of the Au/dielectric inter- 
face [curve (1) in Fig. lb] and the 2-nd band SPP at 
the Au/air interface [curve (3) in Fig. lb] with an energy 
around 1.55 eV. On the other hand, the probe beam inter- 
acts with the SPPs in the 2-nd band of the Au/dielectric 
interface at 1.63 eV around k = 0. Zero order reflection 
and transmission spectra of the investigated structure 
measured using spectrally broad halogen lamp demon- 
strate Wood anomalies 23 related to the SPPs. Fig. lc,d 
show pronounced resonances with an asymmetric Fano- 
like shape 24 . The SPP energies can be determined from 
absorption spectra and are indicated by the arrows. 

The dielectric constant of Au has contributions from 
conduction and bound electrons and can be written as 



xf 



ixf, 



(1) 



where lu p is the plasma frequency, 7 is the electron scat- 
tering rate, and xf an d ixf are * ne rea l an d imaginary 
parts of the interband susceptibility. While xf i s cen ~ 
tered around the electron transition frequencies, xf 1S 



FIG. 2: Time-resolved differential transmission (a) and reflec- 
tivity (b) signals measured for different polarization configu- 
rations of pump and probe beams. The signals are shifted 
vertically for clarity (zero levels are shown by dashed lines). 
The signals for (p,p), (s,s), and (s,p) configurations are mul- 
tiplied by 10. Pump fluence here and for the experimental 
results in all subsequent figures is 0.5mJ/cm 2 . 



spread in frequency and gives noticeable contributions to 
e m even far from the (i-band transition in gold 25 . Op- 
tically induced ultrafast changes in 7, xf \ an d ixf are 
consequences of modifications of the electron and lattice 
properties of gold, caused by electromagnetic energy ab- 
sorption. 

Intraband absorption gives rise to a strongly non- 
equilibrium electron distribution around the Fermi 
energy 16 . Due to electron-electron scattering, electron 
thermalization takes place leading to a hot Fermi dis- 
tribution. Simultaneous to the internal thermalization 
electron-phonon interactions transfer the energy to the 
lattice in a matter of 1 ps. This process leads to the 
transient changes in xf ■ The values of Axf are approx- 
imately proportional to the integrated absorption of the 
laser pulse and evolves in parallel with the electron-lattice 
thermalization 18 . Intraband absorption also changes the 
scattering rate 7. However the temporal behavior of A7 
is different from the one for Axf- While Axf reaches 
its maximum as soon as the pump pulse is finished, A7 
increases during the first hundred femtoseconds and then 
decreases on time scales much larger than 1 ps 17 . This 
is due to the fact that at the earlier stages A7 is mainly 
due to the heating of the electron and phonon systems, 
while for the subsequent stages it is dominated by en- 
ergy relaxation of these systems via electron-phonon and 
phonon-phonon interactions 15 . Along with the excitation 
of the conduction electrons, there is also interband ab- 
sorption as long as ftwi + hw2 > fi^Jib which is satisfied for 
our case. It leads to ultrafast changes in xf persisting 
shorter than a hundred femtoseconds 17 . 

Plasmonic crystals possess extraordinary optical trans- 
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mission and, therefore, transient changes in the gold per- 
mittivity can be monitored not only in reflection but also 
in transmission. Since SPPs are TM-polarizcd they can 
be excited only by p-polarization of the pump beam or 
by s-polarization of the probe beam. The strongest tran- 
sient signals of differential transmission AT/T and re- 
flection AR/R are expected in the case when both pump 
and probe excite SPPs. Indeed, pump SPP allows signif- 
icant electromagnetic energy absorption in gold, whereas 
probe SPP causes Fano resonances in T and R spectra. 
The latter increases the sensitivity of the transient signals 
due to the changes in the optical properties of gold. Ex- 
perimental results prove this statement unambiguously 
(Fig. 2). As expected, the signal is strongest for the fully 
plasmonic configuration (p, s) (p- and s-polarizations of 
pump and probe, respectively). Compared to this con- 
figuration the signal strongly decreases for all other po- 
larization combinations of pump and probe 28 . 

The signal in the totally non-plasmonic case of s- 
polarized pump and p-polarized probe can hardly be 
seen. In the (p,p) case the SPP is excited by the pump 
beam but the changes in reflectivity and transmission are 
more than one order of magnitude smaller. This is ex- 
plained by the weak sensitivity of both transmitted and 
reflected probe beam due to absence of resonances in the 
probed spectrum, and demonstrates the importance of 
SPP for transient optical property sensing. Finally, in 
the (s, s) configuration a plasmonic resonance is present 
in the T and R spectra but the signal is still as small as in 
the (p,p) configuration. This manifests the importance 
of SPP excitation by the pump beam, which provides an 
efficient trapping of the incident photons and their con- 
secutive absorption in the gold. 

The transients AT/T and AR/R versus time t can be 
well described by a sum of two contributions 
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where I — R,T, j = R, T, tj is the relaxation time, and 
(7 = 7£>/\/21n2. Here, the first term corresponds to a 
fast instantaneous component, which follows the auto- 
correlation function of the laser pulses. This component 
is also present in AR/R from flat gold when the pump 
and the probe are incident from the magnetic film side. 
Therefore, we attribute this signal to the non-linear op- 
tical response from the magnetic film related with opti- 
cal transitions in octahedrally coordinated Fe 3+ ions 26 . 
The next contribution is described by the second and the 
third terms in Eq. 2. It manifests a fast increase of the 
signal followed by a monotonous decay to a plateau. This 
contribution is associated with optically induced excita- 
tion of electrons and their energy relaxation in gold. The 
plateau is due to much slower relaxation processes in the 
phonon system. This contribution can be described by 
a multi-exponential decay convoluted with the apparatus 
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FIG. 3: Left panel: Transient transmission (a) and reflection 
(b) at different probe energies ^2 = 1-53 eV, 1.57 eV, 1.61 eV, 
1.68 eV. The solid lines are fits with Eq. 2. Right panel: 
Spectral dependencies of the amplitudes Aj (c) and the decay 
times Tj (d) of the transient signal in transmission (j = T, 
solid symbols) and reflection (j = R, open symbols) evaluated 
from fits of the experimental data. Vertical arrows indicate 
the energy of the probe SPP. 



function. For simplicity, we consider this optical response 
as a single exponential decay with amplitude Aj and de- 
cay time Tj reaching a plateau at the level Bj. In the 
remainder of this work we concentrate on transients re- 
lated to the fast exponential decay. 

The polarization properties (see Fig. 2) directly show 
the major role of SPPs for efficient energy absorption and 
high probe sensitivity. Consequently, the spectral depen- 
dence of the transient signals should also reflect resonant 
behavior and have extrema in the vicinity of the plas- 
monic resonances. In order to check this, we introduced 
an electrically tunable interference filter (20meV band- 
width) in the reflected/transmitted probe beam, directly 
in front of the photodctector. Thereby, we were able to 
probe the evolution of the SPP resonances after absorp- 
tion of the pump pulse (see Fig. 3 a,b). Both AT/T 
and AR/R signals change sign around 1.60-1.66eV. This 
energy region clearly corresponds to the SPP at the 
gold/BilG interface probed in the experiment (sec curve 
(1) in Fig. lb). 

In order to extract information about changes in 7, 
xf, and X2 we analyze the spectral dependence of the 
optical response using Eq. 2. The spectral dependencies 
of the coefficients At and An show resonances with a 
sign change around the plasmonic resonance of the probe 
(Fig. 3 c). Such behavior can be explained by a shift of 
the Fano resonance along with a broadening, which indi- 
cates that both real and imaginary parts of e TO change. 
Surprisingly, the relaxation times tt and tr also vary 
with energy (Fig. 3 d). This indicates that the relative 
contributions of different parts of the transient permit- 
tivity in AT/T and AR/R are energy dependent. 
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FIG. 4: Spectral dependencies of At (a) and An (b) obtained 
by RCWA modeling, to have the best correspondence with the 
experimental data. (c,d) Spectral dependencies of different 
contributions related with A^i 6 (i = 1, solid lines) and Ax}z 
(i — 2, dashed lines) in accordance with Eq. 3. 



In the spectral region, which is addressed in experi- 
ment, A7 influences only the imaginary part of e m . As a 
result, A7 and Ax\ b should have similar spectral depen- 
dencies. However, as discussed above Ax\ b relaxes much 
faster than A7 docs. Since A7 relaxes during times much 
longer than several picoseconds (which is out of the ob- 
servation range) it causes a plateau in the transient sig- 
nals described by Bj in Eq. 2. On the contrary, Ax\ h 
relaxes much faster and contributes to the exponential 
decay with the amplitudes Aj in Eq. 2. Similarly, Axf 
also contributes to Aj. Thus, in the perturbation regime 
the Aj can be written as linear combination of Axl & and 
with frequency dependent coefficients C\ j{uS) and 
C 2 ,»: 



x\ 



Ax? 
xf ' 
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In our experiments the pump energy is far from the d- 
band transitions in gold so x % \ and X\ b can be consid- 
ered constant for the energy range of interest. Spectral 
behavior of the coefficients in eq 3 can be determined 
by modeling the T and R spectra applying a rigorous 
coupled waves analysis (RCWA) 27 . Equation 3 gives 
the best fits of A 3 {oo) for Axf/xt = 0-4 x 1CT 2 , and 
Axf/X2 = 2-6 x lO- 2 (Fig. 4). 

It is seen from Fig. 4 that the relative contributions 
of A\i and A X f determined by Ci j (uj) A X f /x? and 
C2,j{u)Ax 2 b ' lx % 2 are energy dependent. In the energy 
range of 1.45 - 1.60 eV, the term C 1>T (uj)Axf /xf pre- 
vails in absolute value over the term C , 2,t(^)Ax 2 6 '/xf 
and the term C2,r(uj)Ax2 /xf prevails over the term 
Ci tR (u)Ax?/Xi' At higher energies of 1.60 - 1.66 eV 
these relations get reversed. Since different contributions 
have different relaxation times this leads to the energy 



dependence of Tj. As discussed above, Ax\ relaxes in 
less than a hundred femtoseconds, which diminishes Tj. 
This explains the experimentally observed drop of tt for 
larger energies where Ax\ b contributes in Eq. 3 stronger 
(see Fig. 3 (d)). The decrease of tr for smaller ener- 
gies can be explained similarly by stronger contribution 
of A x f. 

Furthermore, the significantly different temporal be- 
haviors of A7 and Ax\ allow us even to discriminate 
their input in the differential signals in spite of the sim- 
ilarity of their spectral contributions. The coefficients 
Bj in Eq. 2 give a A7/7 = 3.0 x 10 -2 corresponding to 
Aj/Axf = 1.15. This nicely agrees with the data in , 
providing A^/Ax? ~ E 

In conclusion, we show that the optical properties of 
plasmonic crystals can be tuned by femtosecond laser 
pulses. This opens a nice possibility for the control of 
the SPP dispersion in the plasmonic crystal. As the re- 
sult SPP propagation through the plasmonic crystal is 
modified and corresponding differential transmission and 
reflection of the order of 3 % is observed. Moreover, the 
response time of plasmonic crystals can be varied widely 
from 200 to 800 fs by tuning the relative spectral posi- 
tions of the probe and the SPP resonance. The results 
are important for operation both in near and far optical 
fields. In the near field it provides the possibility for ul- 
trafast all-optical control of SPPs propagating through a 
plasmonic crystal. In the far field it allows one to modify 
the reflection and transmission spectra of the plasmonic 
crystal, thus making this material of great importance 
for nanophotonic applications. Moreover, the plasmonic 
crystal is also interesting for fundamental research, as it 
provides Wood's anomalies which allow clear distinctions 
among the contributions of the different gold permittiv- 
ity components. An analysis of transient optical spec- 
tra along with a rigorous solution of Maxwell's equations 
provide sufficient information for determining magnitude 
and dynamics of Axi, A%2, and A7. 
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Appendix A: Sample preparation 

The magnetic part of the magnetoplasmonic structure 
is a 2.5 (im thick bismuth-substituted rare-earth iron gar- 
net film of composition Bio.4(YGdSmCa)2.6(FcGcSi)50i2 
grown by liquid phase epitaxy from Bi203 : PbO : B2O3 
melt on the Gd3Ga5 0i2 substrate with orientation 
(111). The film possesses a uniaxial magnetic anisotropy 
in the direction perpendicular to the film plane. The 
specific Faraday rotation is 0.46deg//^m at a wavelength 
of 633 nm. The magnetoplasmonic sample of structure 
shown in Figure la was fabricated by the thermal 
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deposition of the gold layer on the bismuth-substituted 
rare-earth iron garnet film and subsequent electron beam 
lithography combined with the reactive ion etching in 
Ar plasma. The sample was characterized by AFM 
and SEM imaging. The following grating parameters 
were obtained: Gold layer height h = 120 nm, period 
d = 595 nm and the air groove width ?• = 110 nm. 



Appendix B: Experimental technique 

A Ti-Sapphirc self-mode locked laser with a pulse 
repetition rate of about 80 MHz was used as the source 
of ultra-short optical pulses with the center wavelength 



around 800 nm and a bandwidth of 80 nm. The laser 
beam was sent through a pulse shaper and a compressor, 
allowing for an overall time resolution of about 40 fs. 
Subsequently, the laser beam was split into the pump 
and probe beams. The pump beam passed a mechanical 
delay line and was modulated with a mechanical chopper 
at a frequency of 1-5 kHz. Both pump and probe beams 
were focused onto the plasmonic crystal with a reflective 
microscope objective into the spot of about 5 /im 
diameter. The intensities of pump and probe pulses were 
set to 500 and 50 /iJ/cm 2 , respectively. The transmitted 
or reflected probe beam was homodyne detected with a 
balanced photodiode and a lock-in amplifier. 
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